Principles of Cosmology 
Useful (and semi-useful) references: 
Most cosmology textbooks discuss the basic cosmological equa- 
tions for A = 0 cosmologies. For material with A # 0, one can 
look at: 
Peebles, 1993, Principals of Physical Cosmology 


Hogg 2000, “Distance Measures is Comology,” astro-ph/9905116 


Ned Wright’s Cosmology Calculator: 
http://www.astro.ucla.edu/~wright/CosmoCalc.html 


The basic assumptions of Cosmology are 








1) The Universe is homogeneous — every observer sees the 
same thing. 

2) The Universe is isotropic — there is no preferred direction 
in the universe. 


The implication of these two statements is that the universe must 
either be static, or have purely radial motion. (For example, if 
the universe rotated, then the preferred axis of rotation would 
violate the assumption of isotropy.) 


Let's assume that the universe is dynamic, and let 


u= the co-moving coordinates of an object 
R(t) = the motion (expansion or contraction) of the universe 
(= the measured distance to an object 


From these definitions, the observed velocity of a galaxy is 
Pas de d 
dt dt 


Now, neither R nor u is observable, but £ is. So, let's substitute 


¢ for u 
v=R (3) = (È) = H (1.02) 


Here, we have defined the variable H(t) = R/R, with the units of 
inverse time. (The unobservable unit of scale attached to R has 
cancelled out!) H(t) is the Hubble parameter; its value today, 
Ho=Ro / Ro, is the Hubble Constant. H(t) is the fractional rate 
of expansion of the universe. 


(Ru) = Ru (1.01) 





The Newtonian Universe 





Consider a universe where the Newtonian laws of gravity apply. 
Let's pick an arbitrary position in the universe, and call it the 
center. Now consider the motion of a shell of material a distance 
R from this center. (R can be the extreme limit of the universe, 
or, if you wish, it can be just a small region.) 





Recall from freshman physics 
that if the universe is homo- 
geneous and isotropic, mat- 
ter outside the shell will have 
no effect on its motion: the 
deceleration of the shell only 
depends on the matter inte- 
rior to it. 


So the deceleration of the shell is 





If we multiply each side by R and integrate over time, we get 
i GM . 
faRa=- | R2 Rdt (1.04) 


Since R is the derivative of R, and R is the derivative of R, both 
integrals are easy (of the form f udu). The result is 


1., GM 

yA sa e E (1.05) 
where E, the total energy, comes from the constant of integration. 
Note that this is nothing but energy conservation. If the potential 
energy of the universe is greater than its kinetic energy, E < 0, 
and eventually there will be a collapse. If the kinetic energy is 


greater than the potential energy, E > 0. 





Now let’s parameterize how fast the universe is decelerating. 
Again, since we can't measure the size of the universe, we must 
somehow make the units of R disappear. If we take the energy 
conservation equation, multiply it by two, and divide by R”, then 


12 I (1.06) 
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or, since R = —GM/R?, 


E ù 
R 2R 2E 
= ea Í; 
(y ia 
If we multiply and divide the second term in the equation by 
R/R squared, then 


RN" 2R/R\?(R\ 2E 
(È) +25 (5) (È) => (1.08) 
RN BR\ (R\ 2E 


We can now define a dimensionless deceleration parameter 


RR 
q(t) = E (1.10) 
Equation (1.09) is then 
2E 
aT c g= i (1.11) 
or AL 
(1 — 2q) (1.12) 


From this definition, it is clear that if q < 1/2, E > 0 and the 
universe is unbound. If q > 1/2, E < 0, and we have a bound 
universe which will collapse. A value of q = 1/2 implies E = 0, 
which is a critical universe. 


We can also write (1.12) in terms of the density of the universe. 
If we start with the acceleration term 





and substitute density for mass 
M 
be 1.13 
Then ; 
4 
i i (1.14) 
R R R 3 R 
which implies that 
3 
= — H°q 1.15 
aia E (1.15) 
For a critical E = 0 universe, q = 1/2, so in this case 
3 
c= —— H° 1.16 
pom ea (1.16) 
If we let h = H)/100 km s”?, then at the present time, 
= 2 = 1.9 x 107 29h? g cm? (1.17) 





We can parameterize the universe in yet another way. If pe is 
the critical density of the universe, then the cosmological density 
parameter 


Q(t) = p/pe (1.18) 
Obviously, if the universe is bound, Q > 1, but if it’s unbound, 
Q < 1. And just as obviously, Q = 2q. 


The Age of the Universe 


What is the age of the universe as a function of Ho and go? To 
compute this number, we can start with energy conservation 


S ene (1.05) 


and as before, multiply through by 2/R? 


(È) _2GM _2E tw 





R R R2 


The mass of the universe (which is presumed to be constant) can 
be evaluated using the universe’s present size and density, t.e., 


4 
M = 37 Ro Po (1.19) 
So 
R\ 8 Ro\} 2E 
0 
“ E pe = — 1.2 


First, consider an empty (Milne) universe. In this case, po = 0, 
so equation (1.20) becomes 


R\ 2E 
(È) = (1.21) 


R = 2E)? (1.22) 


or 


This is easily integrated from t = 0 to the present to get 


R(t) = 2E) 1 (1.23) 


The age of the universe is therefore 


t= En (1.24) 


This can't be evaluated, but we can substitute R for (2E)!/? 
(using 1.22). So 
R 1 


t = — = — 
R H 


So, for a Milne universe, the present age is to = 1/Ho. 


(1.25) 





Now consider an Einstein-de Sitter (critical) universe. In this 
case, E = 0, and equation (1.20) becomes 


fto fi) o 


If we multiply through by R?, then the equation becomes 


8 
R — 37 GpoR R == 0 (1.27) 
or 
l 8 1/2 
R= EI Roe (1.28) 
The solution of this differential equation has the form 
R = at?” (1.29) 


where 
1/3 


a = {6nrGpo Ri } (1.30) 


(The confirmation of this is left as an exercise to the doodling 
student.) 


Now, to calculate the age of the universe, start with R = at?/3 
and R = (2/3)at-!/3, We can again get rid of the unobservables 
by taking a ratio 


R (2/3)at-!/3 21 
ef a a 1.31 
(È) at2/3 3t Ba 
S 2/R DA 
tS (=== 1.32 
ia 0.82) 
and, at the present time 
Del 
to = -— (1.33) 


Cosmological Redshift 


What effect does the universal expansion have on light? To see 
this, first consider a particle with velocity Y moving past an ob- 
server at point 1 on its way to an observer at point 2. By the 
time it gets there, the universe has expanded; specifically, while 
the particle has traveled v(t)dt, the universe has expanded by 
H v(t)dt. Consequently, the observer will measure the particle’s 
velocity to be 


v(t + dt) = v(t) — (È) v(t)dt (1.34) 
A bit of algebra yields 


n ES 


ps d dR 
U 
_ _dk 1.36 
v R ) 


which simply integrates to v = R7?. The same argument applies 
to the frequency of a photon, yielding 


yt) = (1.37) 


If we now define redshift as 


Ve Abe 


(1.38) 





| 





1 = 
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then 
a= BGA (1.39) 


Einstein’s Universe 


Newtonian gravity gave the energy of the universe as 


(È) tram (EE can 


In the relativisitic case, the equation is very similar 


s 2 
R Ro = ke Ae? 


Here, k is a constant which defines the curvature of the universe. 
Note that k only appears when divided by R? (which is not mea- 
sureable). Thus, to simplify things, we can rescale R (which 
is not measureable) such that k is either 0, 1, or —1. The value 
k = 1 represents a universe with positive curvature, like a sphere; 
k = —1 reflects negative curvature, like a saddle. In the critical 
case, where k = 0, the universe is flat. 


The other new variable A is a Cosmological Constant. You can 
think of it as a pressure term which supplies a new repulsive (or 


attractive) force that is directly proportional to distance, 


R=AR (2.03) 


In a relativistic universe, you must also take into account that 
the distance between two co-moving points is not simply du. In 
normal space, the distance between two points is 


ds? = du? = dx? + dy? + dz? (2.04) 


The total distance, s is then found by integrating along the path. 
In space-time, the interval between these two points is 


ds? = edt? — du? (2.05) 


where the path in space and in time must be integrated. Finally, 
in the case of an expanding universe, the spatial distance between 
two points is Rdu, where R is the size of the universe at the time 
of the measurement. So the distance between two points is 


ds? = códt?* — R° du? (2.06) 


This is called the Robertson-Walker metric. Note that for light 
(which travels at the speed of light), ds = 0. 





An additional complication comes from the fact that space-time 
is not necessarily flat. In Cartesian coordinates 





du? = dz? + dy? + dz? (2.07) 
which, in spherical coordinates is 
du? = dé? + €7d6? + €? sin? 0 do? (2.08) 
But this assumes space is flat. If space is curved, then 
dé? 
Di 2 192 Sr 9 2 
du capa + E“ sin“ 0 do (2.09) 


is the correct expression for distance. Note that if k = 0, the 
equation reduces to normal spherical coordinates. However, if 
k = 1, then space is elliptical, with E < 0, and q > 1/2. Con- 
versely, if k = —1, then we are in the hyperbolic space, with 
E > 0 and q < 1/2. 





Solutions for the Friedmann A = 0 Universe 


The first Einstein Field equation is 


5 2 

R 8 Ro° ke Ne? 

a E ch PE AE 2.02 
12 a ($) RT (2102) 


Now let's take A = 0. If k = 0, then space is Euclidean, and 
(2.02) reduces to 


8 1/2 
R= [raon] Ro (2.10) 


As before, the solution to this equation is the Einstein-de Sitter 
universe 


R = {6Gp R y P 2/3 (2.11) 


If k # 0, however, the solution to the field equation cannot be 
written in closed form. In other words, the solution to 


‘ 8 = 
R? — gr GpoRoR l = —ke? (2.12) 


must be written parametrically. Let us define an intermediate 
variable, 0. In terms of 6, the solution to (2.12) is 


k = 1 (bound) k = —1 (unbound) 
R(0) = a(1 — così) R(0) = a(coshé — 1) 





t(0) = É(0- sino) ig) = (sinh — 0) (2.13) 
C C 
where a 
ae 200 pi (2.14) 


3c2 


Although this parametric form may seem awkward, it is very 
convenient for cosmological calcuations. For example, in a closed 
universe, 


R = a(1 — cos 0) dR = asin 0 de 
t = £(0 — sind) dt = £(1 — così) dé 


So the Hubble parameter, in terms of 0, is 


a/a pre ili cd 


R dt “ (1 così) 
c sind 
ee (2.15) 
a (1 — così) 


Similarly, after a bit of math, one can derive all the following 
relations: 


k= 1 (bound) 
R = a(1 — così) 
dR = asin 0 dé 


a 
t= -(0 — sinó 
S sin 0) 


dt = -0 — così) dé 


c__ sind 

a (1 — cos 9) 
Q _ 1—cos0 
2 


H = 


q = 


sin? @ 
Di) 
Q 


1 — cos 4 


cos 0 = 


1 = ———— 
CARE) 1 — cos@ 
z + cos ĝo 


cos 0 = KETE 


k = —1 (unbound) 
R=a(cosh@ — 1) 
dR = asinh 0 dé 


q= “(sinh 0 — 0) 
C 


= = (cosh 6 — 1) d0 


C sinh @ 
a (cosh 6 — 1)° 
E Q E cosh 0 — 1 
SI a sinh? 9 
2-0) 
bD === 
COS O 
cosh ĝo — 1 
1 SSeS 
nee) cosh 0 — 1 
aie z + cosh ĝo 


(1+ z) 


(2.16) 
(2.17) 


(2.18) 


(2.19) 


(2.20) 


(2.21) 
(2.22) 
(2.23) 


(2.24) 


Cosmological Proper Distance 





There are several ways one can measure distance in astronomy. 
One can compare an observed flux to an emitted flux, and use 
the 1/r? law: such a distance is called a “luminosity distance.” 
Similarly, one can measure the angular size of an object, and 
compare it to its true physical size; this would be an “angular 
size distance.” If the transverse velocity of an object is known, 
then its “proper motion” distance can be found by observing 
its motion on the sky. Finally, and most basically, there is the 
“parallax distance” to an object. 


In cosmology, each of the above distances has a different depen- 
dence on Ho, Qo, and z. Before considering them, however, let's 
calculate the proper distance a ray of light covers in going from 
an object at co-moving coordinate u emitted at time tı, to an 
observer at co-moving coordinate u = 0 at time to. To do this, 
we start with the Robertson-Walker metric 


ds? = cédt? — R°du° 


de? 
1 — kê? 





= códt?* — R° l + £°d0° + €? sin? sas (2.25) 


For light travel, ds = 0, and, since the light path is purely radial, 
d0 = dd = 0. So the Robertson-Walker metric simplifies to 


de? 


ETO > 
c dt =: 





(2.26) 





or, if we take the square root of both sides and move R to the 
left-hand side of the equation, 


d 
Í t= $ 


a TRE (2.27) 





The proper distance between (u,t,) and (0, to) is therefore given 


° T Se E (2.28) 
J, R "i (1 — ké?)1/2 | 


Now consider an Einstein-de Sitter (flat) universe. If k = 0, then 
the right-hand side of the equation is trivial. Furthermore, if we 
use 

R=at?/ (1.29) 


the left-hand side is only slightly harder. Therefore (2.28) reduces 
to 





3c 1/3 3c 1/3 
si Mb = 2.29 
ao a u (2.29) 
Now recall that the definition of redshift (1.39) 
Roya” fi 
1 == = | — 2.30 
( +2) Ri at?’ ti ( ) 
With this substitution, (2.29) becomes 
u= ili a (14237112) (2.31) 
a 
Finally, note that R = (2/3) at71/3, so 
2 
u= = | (1422) (2.32) 
Ro 


and, through the definition of Ho = Ro / Ro, 


= _ 2c = 79 
dy = Rot = 77. {1 (142) \ (2.33) 


Of course, the above equation is only good in the k = 0 case. In 
the general case for a A = 0 universe, 
2c 


dy = ——_—~ 


{oz + (Q — 2) (202 pa 1 \ (2.34) 


Cosmological Luminosity Distance 


Most distances in astronomy are computed via the 1/r? law for 
light. However, for objects at cosmological distances, there are 
additional considerations. The flux of n photons per unit time 
emitted from a source at redshift z can be expressed as 


n-hve photons - energy 








F= 2.35 
dte time ( ) 
Meanwhile, the flux observed is 
n- hvo 1 
= - 2.36 
re (2.36) 


The value d, is the actual distance the photons covers, 1.€e., the 
proper distance. To relate vo, dto to ve, dte, first consider the 
definition of redshift 


Vo = Ve/(1+2) (1.38) 


Next, note that dtp # dt.: because the source is moving, time 
dilation occurs, causing us to measure the source’s clocks to be 


slow, i.e., 


te 
E 2 (2:97) 


{1- way 


In addition, the interval between two pulses will appear longer, 
because the distance between us and the source is ever increasing. 
(A second pulse has a longer distance to travel.) This extra time 
19 
ie A (2.38) 
{1- (w/e) ) 


Thus, 


MATTEI 2 1 — (v/c) 


or 
dto = (14+ z)dt. 


So, with these two additional terms, we have 


nhvo 1 nhv. 1 1 F 


I=" 27 +222 (1422 





1+(0/c) i ca LN =(1+2) 


(2.39) 


(2.40) 


(2.41) 


The luminosity distance is therefore related to the proper dis- 


tance by 
di = dy (1 + z) 


(2.42) 


Angular Diameter Distance 





Consider a standard galaxy with lin- 
ear size r at redshift z. Under normal 
Euclidean geometry, the angular size 
the galaxy subtends would be inversely 
proportional to distance. However, in 
an expanding relativistic universe, the 
calculation is a bit more complicated. 


Again, let's begin with the Robertson- 
Walker metric 


V 


ds? = c°dt° — R? du? 


de? 
1 — kë? 





= ¢c’di* — R? l + €7d6? + €? sin” gag? | (2.43) 


If the galaxy is in the plane of the sky, then the radial distance 
to both sides of the galaxy is the same. Thus, dé = 0. Similarly, 
since both sides of the galaxy are being observed at the same 
time, dt = 0, and the separation of the two sides of the galaxy is 
simply fds = s. Finally, we will choose our coordinate system 
such that the angle subtended by the galaxy is entirely in the 0 
direction, so that d@ = 0. Thus 


ds? = —R*du? = — R°&°d0° (2.44) 
or, if we choose the coordinates to make @ positive 


0= + (2.45) 


Let's get rid of the R (the size of the universe at the redshift of 
the galaxy) and substitute Ro (the size of the universe today), 
using 


Then Ass) 
S 8 E 
d= a ne" (2.46) 


Now recall that Ro€ = Rou is the proper distance. So 


_ 8(1+2) 
= th (2.47) 


In other words, the angular size distance 
da =dp(1+2)1=dz(1+ 2) (2.48) 


Note what this means for surface brightness. From dimensional 
analysis, the surface brightness of a galaxy is given by 


ES 


= 23 


(2.49) 


where f is the observed flux, and 0 the observed angular size. 
Since f is related to the intrinsic flux, F, via the equation for 
luminosity distance, and @ is related to the object’s true size, r, 
via the angular size distance, then 


FO Fda F ni 


Objects undergo a rapid decrease in surface brightness at high 
redshift! 


The Cosmological Volume Element 


For many calculations, you need to know how the volume of the 
universe changes with redshift. In general, the volume element 
should have the form 


dV = r° sin 0 dr d0 do (2.51) 


dV 


The angular variables 0 and @ are unaffected by cosmology and 
act the same as they would in Euclidean space. Since the volume 
shell is concentric with the earth, the radius vector, r, is simply 
the angular size distance. However, r is not a very useful variable 
— the observable variable is z. Thus, a more useful way to write 
the volume element is 


d 
dV =r sin 97 dzd0dé (2.52) 


Note that here, because r is changing with z, use of the angular 
diameter distance is not appropriate for the derivative; for this. 
Instead, we need the proper distance at the distance of the shell, 
i.e., Rdu. 


To calculate the volume element, we start with the Robertson- 


Walker metric 
ds? = cdt? — R° du? (2.06) 





which, for light, means Rdu = cdt. Since we are interested in 
how Rdu changes with z, we can translate dt to dR using 


a e oar (ss 


and convert dR to dz via 





Ro Ro 
== = dk=—--——d 2.54 
R (1+ z) R +23 Z (2.54) 
This gives us 
C C Ro C 
= db dz = —_-d 2. 
Ne Rae aaa 9) 


Now we must relate the Hubble parameter, H, to today’s Hubble 
Constant. Recall that in the case of a closed universe, 


c sin 0 
H < 2.20 
a (1 — cos 0)? E 
SO RE RS 
H _ sin (1 — cos ĝo) (2.56) 


Ho  sin0o (1 —cos0)? 


From (2.23), the second term on the right-hand side is simply 
(1+ z)?, while the first term can be found through (2.24) to be 


sinó — (Qoz + DE 


S Gee ee 


This gives 


H 
= = (Noz + 1)? (+2) (2.58) 
U 
So 4 
Rdu = om i = 53 


H(1 +2) Ho (Qoz + 1)1P (1 + 2)? 
The volume element is therefore 


C 


dV = r4 
Ho (Qoz +1)? (1 + 2)? 


sin 0 dz dô do (2.59) 


Cosmology with a Cosmological Constant 


If there is a cosmological constant (or vacuum energy), the basic 
equations for distance and lookback time are somewhat changed. 
First let’s look at the Einstein Field equation 


. 2 
R 8 ke? Ac? 
AN o a a 2.02 
12 STD RS + (2.02) 


In order to have a flat universe with no cosmological constant, 
the critical density of matter must be 


3772 


= — 1.1 
STG di 


Pc 


If A is non-zero, however, the density needed for a flat universe 
is decreased 
3H? = Ac? 


.01 
87G GOL) 


Pp Z 


Because the “critical density” is no longer a well-defined quantity, 
what is commonly done is to keep (1.16) and define three new 
quantities: 





8TG 

Opie a 02 

Ac? 

On. E l 

\ 3H (Gila 

kc? 
Qr = = .04 
k R? H2 (9:05) 

with 

Qu + Ma tO, =1 (3.05) 


Note the meaning of each term. Qm is the dimensionless matter 
density of the universe, (2, is a dimensionless quantity repre- 
senting the “energy density” of the cosmological constant, and 
Qg is a dimensionless false density that describes the universe’s 
departure from flatness. In an inflationary universe, Q% = 0, and 
Qu + (nr = 1. 


Using these variables, the equation for the expansion of the uni- 
verse can be re-written and simplified. 


A 2 
R\ | RON” ke? A 


Multiplying through by 1/H@ yields 


; 2 
RoN { R SARE 0 ke? Ro i A? 
_— = TT — —s 
Ro) Rl 3R L R RH?) (R) 3H 


(3.06) 





or, after dividing through by (R0/R)?, 


PR (Y er 


Finally, if we substitute for Q, using (3.05), and let x = R/Ro, 
we get 


Cal (de) ie € z 1) +9 (221) (3.08) 


This equation describes the universe as a function of three univer- 
sal constants, Ho, Qm, and Qa. From this (with a little algebra), 
we can also derive qo 


1 


Lookback time with A 
Using (3.08), It is fairly straightforward to calculate lookback 


time with A. Since x = R/Ro = (1+ 2), (3.08) can be re- 
written as 


d(1+2) S ò 2° +22 
—————_ | |=H<{1+Qmz-Qx-——T_ 
( dt ) D CET 


(1+ ye = Ho(1+2) {(1 +2)? (1+ wsi — Qaz(2+ 23 PP 


di 
(3.10) 


This gives 


At = HT Paro? {(1+2)°(1+Qmz)- Qaz(2+ ay? dz 


(3.11) 
Unfortunately, unless QA = 0 or Qm + Qa = 1, this integral 
must be evaluated numerically. Note however, that as Q4 be- 
comes larger, the denominator of (3.11) becomes smaller, and 
the lookback time increases. The additional time reconciles the 
expansion age of the universe with the age of its stars. 


To a few percent, the total age of a universe with A can be 
approximated by 


2 ih O 
A e (3.12) 


ta S 
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where 
Qa = An — 0.301 + 0.3 (3.13) 


(Note that if Q, > 1, you need to replace sinh with sin and flip 
the signs. Note also, that if Qa = 1, the equation is exact.) 


Cosmological Distances with A 


The evaluation of cosmological distances in a universe with non- 
zero A follows closely that of the Friedmann case. For light fol- 
lowing a radial path, the Robertson-Walker metric gives 

C dE 

—dt = —— —; 2.20 

B ui 


If we use (3.10) and manipulate the left-hand side of this equa- 
tion, we get 


ti ti 
Pla] oes 
to R to R Ro dz dt 


T a 
o Ro R dz 
= / = {(1 +2)? (1+Qmz) — Qn2(2 + ay? dz 
o HoRo 
Le 42 tA (Oey) PAE, ee 
Ie 
0 
(3.14) 
Meanwhile, the right-hand side of (2.27) gives 
6 d 
/ a = sinh * U (3.15) 
0 — 


for Qy < 1 (and sin! u for Q > 1). Thus 
d, = Rou 
= Rosinh Ti n {(1+2)? (1+Q2) -Ma2(2+ 2) 4 az} 
i (3.16) 


So if we define 


E= f {(1+2)?(1+Qyz) — Qaz(2+ 2)} TE da 


then 


dp 


= |? sinh {|O,|!/2E} 
0 


€ 
"H, 


0) 
Ho 


SLA sin (QUE) 


(3.17) 


Qg <0 
Oj =0 


Or > 0 
(3.18) 


